Enhanced non-quark-antiquark and non-glueball Nc behavior of hght scalar mesons. 
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We show that the latest and very precise dispersive data analyses require a large and very unnat- 
ural fine-tuning of the 1/Nc expansion at A^c = 3 if the /o(600) and -ft'(800) light scalar mesons are 
to be considered predominantly qq states, which is not needed for light vector mesons. For this, we 
use scattering observables whose 1/Nc corrections are suppressed further than one power of 1/Nc 
for qq or glueball states, thus enhancing contributions of other nature. This is achieved without 
using unitarized ChPT, but if it is used we can also show that it is not just that the coefficients of 
the 1/Nc expansion are unnatural, but that the expansion itself does not even follow the expected 
1/Nc scaling of a glueball or a qq meson. 
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Light scalar resonances play a relevant role for several 
fields of Physics: For the nucleon-nucleon interaction, be- 
cause they are largely responsible for the attractive part 
[l| (with cosniological and anthropic implications), for 
the QCD non-abelian nature, because some of these res- 
onances have the quantum numbers of the lightest glue- 
ball, also common to the vacuum and hence of relevance 
for the spontaneous chiral symmetry breaking. Moreover, 
they are also of interest for the saturation [2[ of the low 
energy constants of Chiral Perturbation Theory (ChPT) 
Q. However, the precise properties of these mesons, 
their nature, spectroscopic classification, and even their 
existence — as for the 7^(800) or k — are still the object 
of an intense debate. In particular, different models Q 
suggest that they may not be ordinary quark-antiquark 
mesons, but tetraquarks, meson molecules, glueballs, or a 
complicated mixture of all these. The problem, of course, 
is that we do not know how to solve QCD at low energies. 

However, since the QCD 1/Nc expansion is applica- 
ble at all energies, and the mass and width Nc depen- 
dence of qq mesons and glueballs is well known [5|j the 
Nc scaling of resonances becomes a powerful tool to clas- 
sify them and understand their nature. In (y, |7[ some 
of us studied the mass and width behavior of light res- 
onances using ChPT — which is the QCD low energy ef- 
fective Lagrangian — and unitarization with a dispersion 
relation. It was found that the poles of the p(770) and 
K*{982) vectors behave predominantly as expected for 
qq states whereas those jaf the /o(600), also called a, and 
K{800) scalars do not Gj. Still, a possible subdominant 
qq cornponent for the /o(600) may arise naturally at two 
loops f?| within ChPT (less so at one loop), but with a 
mass around 1 GeV or more. 

Of course, all these conclusions rely on unitarized 
ChPT and the assumption that corrections, suppressed 
just by 1/Nc, are of natural size. Since Nc^3 in real life, 
this may not seem as a large suppression, even more when 
the meaning of "natural size" may not be clear for dimen- 
sional parameters. For that reason, unitarized ChPT was 
useful to change Nc, and reveal the 1/Nc scaling, no mat- 
ter how unnatural the coefficients may appear. 

Here we will provide adimensional observables with 
corrections suppressed further than 1/Nc, that can also 



be applied directly to real data at iVc = 3, without the 
need to extrapolate to larger Nc using unitarized ChPT. 

In particular, resonances appearing in elastic two-body 
scattering are commonly identified by three criteria. The 
Nc behavior of one of these criteria — the associatedpole 
in the unphysical sheet — was already studied in [q, Q. 
A second possibility is to define the mass as the energy 
where the phase shift reaches 7r/2, which both for tttt or 
ttK scattering occurs relatively far from the /o(600) and 
K{800) pole positions. This criterion was studied in [g] 
for the /o(600) with a relatively inconclusive result about 
its assumed qq behavior. A more reliable parametrization 
and better data were called for and we will provide them 
here together with more conclusive results. Third, the 
phase increases very fast in the resonance region and the 
mass can be identified with the maximum of the phase 
derivative. All three criteria roughly coincide for nar- 
row resonances, but the most physical definition is the 
latest, since it identifies the resonance as a metastable 
state whose lifetime is the inverse of the width. Note 
that this is the less evident feature both for the /o(600) 
and ^(800) and thus the phase derivative will become 
our preferred observable to test their Nc dependence. 

Let us then recall that partial waves generically scale 
as 1/Nc, except at the resonance mass mji. Actually, 
it has been found Q that if a resonance pole at sr = 
TO^ — imnTji behaves as a qq 5], i.e. mji ~ 0(1) and 
Tn ~ 0{1/Nc), then the phase shift satisfies |23 |: 
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where t{s) is the partial wave, a — Ik/ yfs and fc is the 
meson center of mass momentum. Derivatives are taken 
with respect to s. The 1/Nc counting of the different 
terms in the equations above comes from the following 
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TABLE I: Normalized coefficients of the l/N^ expansion for 
different resonances. For qq resonances, all them are expected 
to be of order one or less. 
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In brief, the corrections in Eqs. ([J) to ([4]) are suppressed 
by a further l/N"^ power due to an expansion on the 
imaginary part of the pole, which scales hke F ^ 1/iVj,. 
As nicely shown in [8[, by expanding separately the real 
and imaginary parts of t""'^, only the 1/7V^"+-'^ powers are 
kept on each expansion, leading to Eqs. ^ and ^. 

Since we are interested in adimensional observables 
whose corrections are suppressed further than just l/Nc, 
we can recast Eqs.(IT]) and ([2]) as follows: 
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Note that we have normalized each equation and ex- 
tracted the leading 1/Nc dependence so that the coef- 
ficients a and b should naturally be 0(1) or less. It is 
relatively simple to make a and b much smaller than one 
with cancellations with natural higher order 1/Nc contri- 
butions, but very unnatural to make them much larger. 

Now, in Table U we show the resulting a and b for the 
lightest resonances found in tttt and ttK elastic scattering. 
Before describing in detail the calculations, let us observe 
that for the p(770) and K*{892) vector resonances all 
parameters are of order one or less, as expected for qq 
states. In contrast, for the /o(600) and ^(800) scalar 
resonances we find that all parameters are larger, by two 
orders of magnitude, than expected for qq states. This 
is one of the main results of this work and make the qq 
interpretation of both scalars extremely unnatural. 

Let us now describe in detail our calculations and their 
different degree of precision and reliability. As com- 
mented above, the /o(600) "Breit-Wigner" mass was al- 
ready studied JS] using Eq. ([l}, but no conclusion was 
reached on whether the deviations were consistent with 
the 1/Nc suppression or not. This was partly attributed 
to the limited reliability of the conformal parametrization 
or unitarized ChPT — whose phase never reaches 7r/2- 
used in [8|. To overcome this caveat we are now using 
the recent, very precise and reliable output of the data 
analysis in \§\ constrained to satisfy once subtracted cou- 
pled dispersion relations — or GKPY equations — as well 



as Roy equations, which is therefore model independent 
and specially suited to obtain the /o(600) pole [l3]. Note 
that this analysis incorporates the very recent and re- 
hable data on K^ decays from NA48/2 [l^, UII, which 
is a key factor in attaining high levels of precision. The 
analysis in [9] is also in good agreement with previous dis- 
persive results based on standard Roy equations 12]. We 
have followed the same rigorous approach for the p(770), 
although, being so narrow, the conformal unconstrained 
data analysis and the lAM yield very similar results. The 
uncertainties we quote for both the /o(600) and p(770) 
cover the uncertainties in the output of the dispersive 
representation. 

In this work we also deal with strange resonances in nK 
scattering. For the scalar K{800) we have also used a rig- 
orous dispersive calculation, namely, that in [13|, which 
uses Roy-Steiner equations to determine the isospin 1/2 
scalar channel of ttK scattering, although this time we 
can only provide a central value. Note, however, that 
the value of m^ obtained in that analysis is located be- 
low threshold, so that the phase shift is ill defined at 
ttir. Nevertheless, we have been using the mn mass def- 
inition to allow for an easier comparison with [8|, but 
the definition ^/sr^ = m — iV /2 is equally valid and is 
actually the standard choice used in the context of scalar 
mesons. Moreover, the Nc scaling of Eqs. ([!]) and ([2]) 
does not change if we evaluate the quantities at s = m^, 
instead of rn^R, since m^ differs from m^ in F^/4, which 
is OiN-"^). Thus, the values for the X(800) in Table U 
correspond to this choice. For the vector K*{%92) there 
are no very precise purely dispersive descriptions of the 
existing data and we therefore rely on a single partial 
wave dispersion relation and SU(3) ChPT to one-loop 
to determine its subtraction constants (this is known as 
ChPT unitarized with the single channel Inverse Ampli- 
tude Method (I AM) [Ijl), which we will briefiy explain 
in the next section. We have applied the same method to 
the p(770) and the results lie within 50% of their central 
value when using the GKPY dispersive representation. 
Since the X*(892) is narrower than the p(770), the lAM 
is likely to provide a better approximation than in the 
/3(770) case, but even with that 50% uncertainty, it is 
enough to check that the a and b parameters are smaller 
than one. 

There is, of course, another way of interpreting our 

results, which is that due to the large 1/Nc coefficients 

of the /o(600) the series simply does not converge. In 

xirticular, Eq.(IT|), which was thoroughly considered in 

, is obtained as an expansion of arctan(a;) = x — x^/i... 
In this way we could explain why the a = —0.06 ± 0.01 
coefficient is so small for the p(770): it is simply the effect 
of calculating a = a^/3 with 5 — 0.56lg'Q4, which is now 
naturally of 0{1). We could try the same procedure for 
the /o(600), assuming its series expansion is that of a qq, 
to find a — 9.1, still rather unnatural, but of course, this 
value makes no sense since the whole series would not be 
converging and terms higher than 1/A^^ would become 
dominant. 



This is one of the reasons why despite being only sup- 
pressed by l/A^c instead of 1/N^, we also provide the 
expansion in Eq. (j6|) obtained from the derivative of the 
amplitude. In this case the b/N^ term is not the square 
of a natural l/Nc quantity, i.e., 
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Despite containing a cancellation between two 1/Nc 
terms, its value for the p(770) is rather natural. How- 
ever, once again, the value for the scalars is almost two 
orders of magnitude larger than expected. 

In the previous analysis it is very relevant that the 
width of the resonance is suppressed with additional 1/Nc 
powers with respect to the mass. Actually, it is rather 
straightforward to extend the formalism to study the as- 
sumption that the /o(600) could be predominantly a glue- 
ball, since then niR - 0(1) and Tr ~ 0{1/N^) ^^. 
As a consequence, for the glueball case, the scaling of 
Eqs. (|3]) and dU changes and so does that of ^(m^) and 
S'iml): 
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Much as it was done in Eqs. ([5]) and (O, in order to make 
explicit this further N^ suppression we can define some 
new parameters a' and 6' that should be of 0(1) if the 
resonance was a glueball: 
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Following the same procedure as before we obtain for the 
/o(600), a' = -6800l!^2oo and b' ^ 2080t65o- In other 
words, a very dominant or pure glueball nature for the 
/o(600) is very disfavored by the 1/Nc expansion, even 
more than the qq interpretation. This is because it would 
require even more unnatural coefficients, this time too 
large by three to four orders of magnitude. 

Of course, as we did for the qq case, we could worry 
about the fact that, due to the arctan(a;) = cc — x^/3 -I- ... 
expansion, the a' should have been interpreted as a' — 
a'^/3. But even with that interpretation we would still 
find a' = 27^7, again rather unnatural. Once more, and 
as it happened in the qq case, the b' parameter does not 
correspond to the fourth power of any natural quantity, 
so that its value is genuinely unnatural, disfavoring the 
glueball interpretation. 

Let us remark that in the case of tetraquarks or 
molecules, the width is not expected to be suppressed 
with additional 1/iVp powers with respect to the mass of 
the resonance |15l.ll6l|. Thus, our previous formalism does 



not apply. Furthermore, it is most likely that scalars are 
a mixture of different components. Therefore our results, 
while showing that neither the qq or a glueball are favored 
as dominant components of light scalars, do not exclude 
that these structures could be mixed with other compo- 
nents that would dominate the 1/A^c expansion with a 
different Nc behavior [23 |. 

In summary, we have just shown that if, for the light 
scalar mesons, we study qq or glueball 1/Nc expansions 
as those in Eqs. ([5]), ([6]) and (TTOt . their coefficients come 
out very unnatural, suggesting that these resonances can- 
not be described as predominantly made of a quark and 
an antiquark or a glueball. Note that, contrary to our 
previous works [y, 0], this conclusion has been reached 
from dispersive analyses of data, without extrapolating 
to TVc 7^ 3 using unitarized ChPT. 

However, unitarized ChPT will be used next to calcu- 
late the Ai — 1 observables, in order to show that, for 
scalars, what really happens is that they do not even fol- 
low the 1/Nc expansion of qq or glueball states given in 
Eqs. di]), dH) and pH]) . thus explaining the need for un- 
natural coefficients if a qq or glueball-like expansion is 
assumed. 

The Inverse Amplitude Method: The elastic lAM 
[14| uses ChPT to evaluate the subtraction constants and 
the left cut of a dispersion relation for the inverse of the 
partial wave. The elastic right cut is exact, since the elas- 
tic unitarity condition Imt = cf\t\'^ , fixes lTa.t~^ = —a. 
Note that the lAM is derived only from elastic unitar- 
ity, analyticity in the form of a dispersion relation, and 
ChPT, which is only used at low energies. It satisfies 
exact elastic unitarity and reproduces meson-meson scat- 
tering data up to energies ~ 1 GeV. It can be analytically 
continued into the second Riemann sheet where poles as- 
sociated to resonances are found. In particular, we find 
the p(770) and /o(600), as well as the if* (892) and the 
K{%OQ) resonances as poles in tttt and iiK scattering am- 
plitudes, respectively. 

The dependence on the QCD number of colors, iVc, 
is implemented [Q, u\ through the leading Nc scaling of 
the ChPT low e nerg y constants (LECs), and is model 
independent 0,0, Ezf- Fortunately, for Eqs.® and © to 
hold, only the leading 1/A^c behavior is needed. Note also 
that the lAM does not have any other parameters where 
uncontrolled Nc dependence could hide — as it happens 
in other unitarization methods — so that the lAM allows 
us to check the scaling of the A^ — 1 in Eqs.© and ^. 

The SU{2) I AM: Only the non-strange /o(600) and 
p(770) resonances can be checked, but we can do it uni- 
tarizing with the lAM the corresponding partial waves 
either to one or two loops. We simply scale /at^ -^ 
f\/Nc/'i, the one loop constants as II j^^ — >■ l^Nc/S and 
the two loop ones as ri^N^ ~^ ri{Nc/3y. 

Thus, in the two first columns of Fig.[T]we show, for the 
p{770) and /o(600) resonances, the scaling of the A^ — 1 
both to one loop (upper panels) and two loops (lower 
panels). Note we have normalized them to their Nc — 3 
value, in order to cancel the leading part of the a and b 
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FIG. 1: 1/A'^c scaling of the A; — 1 observables normalized to their Nc = 3 value for light scalar and vector mesons, using 
unitarized ChPT within SU{2) or SU{3) and to one or two loops — 0(p'*) and 0{p^), respectively. 



coefficients and thus extract the leading l/N^: behavior 
of Eqs. (O and ^. For the one-loop calculations we use 
the set of LECs in [18], whereas for the two- loop calcu- 
lation we use the fit D from [l^ [l^ . We have checked 
that similar results are obtained when using other sets of 
LECs in these references or the estimates from resonance 
saturation [2]. 

We can observe that the scaling for the p(770) observ- 
ables overlaps with the expectation for the leading behav- 
ior of qq states. However, in the case of the /o(600) the 
scaling is completely different. To one loop the /o(600) 
observables grow instead of decreasing. Let us note how- 
ever, that, for Nc larger than ^ 10, the /o(600) pole lies 
on the third quadrant of the complex plane. Before that 
happens, the value of m|j becomes less than 4m^ and the 
phase shift has no physical meaning so that Eqs. ([SJ and 
(ini) do not hold. This behavior does not occur to two 
loops. Actually we find again the /o(600) behavior al- 
ready observed in Q, where, for Nc close to 3, the width 
grows as in the one-loop case (and so do the observables 
here), but for larger Nc the /o(600) starts behaving more 
as a qq. Note that this qq behavior appears at a mass 
somewhat bigger than 1 GeV. This was a hint of the 
/o(600) being a mixture of a predominantly non-qq com- 
ponent and, at least, a subdominant qq component with 
a mass much heavier that the physical one, which is the 
one that survives at large Nc . In terms of the A^ — 1 ob- 
servables defined here, this translates into a growth close 
to A^c = 3 and a decrease at larger Nc. Therefore, it is 
not only that the a and b coefficients of the /o (600) are 
too large as shown in the previous section, but that the 
scaling itself does not correspond to a qq state (and even 
less so to a glueball). To two- loops, the p{770) does not 



follow exactly the leading behavior of qq states but de- 
creases slightly faster, which can be naturally explained 
due to subleading effects or to a possible small pion cloud 
contribution. 

The SU{3) I AM: Now we can study the scaling of A^ - 
1 not only for the p(770) and /o(600), but also for the 
ii'*(892) and if (800) resonances, although in this case 
the elastic unitarized amplitudes are available only to 
one loop [2 a[2ll We have now eight LECs, called Li{p), 
that scale 0, [i3l as Li^^^ — >• Li{Nc/5) for i — 2,3,5,8, 
while 2Li — L2, L^, Lq and L7 do not change with Nc- 

In the third and fourth columns of Fig. [T]we show the 
results found using the set of LECs called Fit II in [23] ■ 
Similar results are obtained with Fit I or the estimates 
from resonance saturation in [2] . In the upper panels we 
simply reobtain within the SU{3) formalism the same 
results we obtained for the p{770) and /o(600) within 
the SU{2) formalism to one loop. In the lower panels we 
show the results for the light vector AT* (892), following 
nicely the qq expectations, as well as the results for the 
scalar A'(800), which has a very similar behavior to the 
/o(600), at odds with a dominant qq or glueball nature. 

Summary: In this work we have studied the 1/A^c 
expansion of the meson-meson scattering phase-shifts 
around the pole mass of a qq or glueball resonance. In 
particular we have defined observables whose corrections 
are suppressed further than just one power of Nc, pay- 
ing particular attention to the derivative of the phase, 
which provides a physical and intuitive definition of a 
resonance. Using recent and very precise dispersive data 
analyses we have shown that if we assume a qq or glueball 
behavior for the /o(600) and ii'(800), the coefficients of 
the expansion of such observables turn out unnaturally 



large. This is shown without using ChPT or extrapo- 
lating beyond N,. = 3. Moreover, when using unitarized 
ChPT, we have shown that it is the very 1/Nc scaling of 
the observables which does not follow the pattern of the 
1/Nc expansion expected for qq or glueball states. 
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